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ABSTRACT 

A Majorana type of the neutrino mass matrix induces a class of lepton 
number violating processes. Cross sections of these reactions are given in 
terms of the neutrino mass matrix element, and a semi-realistic event rate is 
estimated. These processes provide mass and mixing parameters not directly 
accessible by the neutrino oscillation experiments. If these processes are 
discovered with a larger rate than given here, it would imply a new physics 
of the lepton number violation not directly related to the Majorana neutrino 
mass, such as R-parity violating operators in SUSY models. 
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Neutrino oscillation observed in SuperKamiokande [[], SNO 0, and Kam- 
LAND 12] experiments have opened a new window beyond physics of the 
standard model. The immediate critical question is to determine the nature 
of neutrino masses; whether they are of Dirac or of Majorana type. Unfor- 
tunately, the data on the neutrino oscillations can not discriminate between 
the types of neutrino masses, as we will see below. In order to get a conclu- 
sive argument on this issue, the investigation into lepton number violating 
processes, inevitable consequence of the Majorana nature, is indispensable. 

If the Majorana mass is verified, it opens up the possibility of generat- 
ing the baryon asymmetry of the universe via leptogenesis scenario H|. A 
conventional lepton number violating process and the one most extensively 
discussed towards this goal is the neutrinoless double beta decay [5^ . Despite 
several ingenious experimental proposals for improved detection of this de- 
cay we believe that some alternative methods to measure the Majorana type 
of the neutrino masse matrix are both useful and very important. In the 
present work we examine a variety of lepton number violating processes for 
this purpose. 

We systematically examine a class of lepton number violating effective 
operators below the Fermi energy scale 1 / ^/Gp, 

II qq qq , (1) 

where / is the lepton doublet and q is the quark doublet having the quantum 
number of the standard model. Existence of this class of operators requires 
a new physics beyond the standard theory, but we do not need to specify 
the new physics. The Feynman diagram that generates this class of effective 
operators is depicted in Figure 1; the cross in the figure is the Majorana 
neutrino mass matrix and the exchanged particle is a weak boson, either 
or Z. When both of / are the electron and q is either a u- or d-type 
quark, the operator gives the neutrinoless double beta decay. We extend this 



2 




Figure 1: The Feynman diagram that generates effective lepton number vi- 
olating operators 



double beta process to all combination of charged leptons of I, such that the 
full neutrino mass matrix element may be explored. 

The strength of these dimension 9 operators at low energy scale of a/s < 
250GeV is of order 



3/2 



10 



-22 



(2) 



leV 'lOOMeV' 

This is an extremely small number compared to the usual weak interaction 
of order Gps, but the nature of neutrino mass may well be examined only 
by these operators. The important point is that once the Majorana nature 
of the neutrino mass is assumed, there is no arbitrary freedom left, both for 
existence and its strength of the class of lepton number violating processes 
we consider below. 

Cross sections of processes discussed below are determined by using the 
matrix element of the Majorana neutrino mass. 
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where [k = 1,2,3) and = e,fi,T) and m^'s are the mass eigenvalues 

(real and positive by definition) . The reaction rate is thus independent of the 
origin of the neutrino mass matrix such as the seesaw mechanism jH]. Note 
that the above combination of neutrino masses and the mixing parameters is 
different from that measured in the neutrino oscillation experiment, in which 
Uak and its conjugate [/^^ appears, thus the Majorana phases, characteristic 
of the Majorana nature, just cancel out in neutrino oscillation. Since the 
Majorana phase factor plays an important role in leptogenesis calculation, 
in general, the investigation into the lepton number violating reactions is 
relevant for the leptogenesis. 

If one of these processes is discovered with a larger rate than given below, 
it means existence of a new class of diagrams not involving the neutrino mass 
matrix and may provide a new feature to the lepton sector such as R-parity 
violating interactions in SUSY models. In that sense lepton number violating 
processes complementary to the neutrinoless double beta decay may provide 
the unique window to determine the mechanism of how the lepton number 
violation occurs. 

With the lepton fiavor mixing, the effective operators above give a variety 
of lepton number violating processes. We consider the following reactions 
which we think relatively easy to explore experimentally; 

(e/i); e- + ^ fi+ + A^-\ (4) 

(/ie); /i- + ^ e+ + (5) 

(e/i); p + A^^fi+ + e+ + iA+lf-\ (6) 

(ee); e~ + e~ {atomic) —^n~ + n~, (7) 

(cyu); + e~ (atomic) + n'^ . (8) 

Here (ajS) means that the process may explore the matrix element ma[3- 
There are corresponding inclusive processes such as e~ + A^ — > /i"*" + X , 
where X is any hadronic state. When A is a light nucleus, it may break up 
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due to a low binding energy, for instance, 

e- + He^ ^ 11-^ + An . (9) 

To reduce the background at low energies to a respectable level, it is 
important to search for unambiguous signatures , for instance without decay 
muons below the vr production threshold for the first process of Q. On the 
other hand, at high enough energies pions may not decay within the detector 
volume, in which case one does not worry about /i^ from the pion decay. It is 
also important to use a high intensity beam and a high density target in order 
to overcome the low reaction rate, for instance the atomic electron target in 
the processes of ((Tj) and (jH)). We also consider muons captured by nucleus in 
© to compensate for the low flux. The other processes we considered and 
not listed here are difficult due to various reasons; the background problem, 
the insufficient beam energy available or planned now, a finite short lifetime 
and so on. 

We consider important processes in turn, and the conversion processes, 
e~ + He^ — > /i"*" + 4n and e~ + — > yU+ + A^~'^ are the first. A related 
conversion process, ^~ + + ^(^-2) ^ ^- _^ _^ ^+ _^ ^{^-2) ^ j^g^g 

been considered in [5] and [7|. The experiment could be done effectively below 
the TT-production threshold rejecting the background from the 7r+ decay, 
hence for this case the incident electron energy should be chosen properly. 
For the 2nd process of nuclear target, the cross section at low energies is 
given by 

- = E--/ = Er^..! , (10) 

where t^^ is the lepton wave function, Vt^j_i, is the nuclear matrix element and 
Rn is an effective nuclear radius, as defined in jHl. Since the final nucleus 
can take any final state, we sum up the final states. Therefore, we use 
the average momentum and energy for yU"*". We neglect the momentum of 
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in t'^'^ and use the approximation to evaluate the nuclear matrix element 
0. We find If^^f.ul^ ^{9v + ^8g\ + 9g\)Z{Z - 1) for 5 = and 
{gy + 2g\ + 9g\)Z{Z — 1) for 5=1. Here, S is the spin of two protons and 
Z is the atomic number. From this, we find with Z = 50 and A = 100 

,.5xlO-c,n^(J^f . (11) 

for the average /i momentum of SOMeV. It may be worthwhile to comment 
that there is no reduction due to the nuclear matrix element in comparison 
with the double beta decay and also there is enhancement due to many 
possible combinations of the proton target inside the nucleus. 

The helium case can roughly be estimated by just taking Z = 2, and we 
find 

In order to compensate for the small reaction rate at low energies, it might 
be useful to go to higher energies despite larger backgrounds. Energy de- 
pendence of the cross section goes like oc with ^/s the CMS energy and 
it is 0[l]Gj^\mefj,\'^ s'^/16tt'^ . A simple parton model neglecting transverse 
momentum may be used in estimating 0[1] factor here. 

The advantage of muon capture listed in (0) is two fold; it automatically 
gives a self-focusing mechanism of the incident particle, capture into the area 
of nuclear size, and the same muon may repeatedly be used as in the case 
of a high luminosity accumulator ring. Thus, a naive computation yields 
a very large event rate, for instance even with a pulsed muon fiux of order 
10^^/ sec much smaller than the highest achieved electron fiux. The problem 
however is that at the same time the first order weak process, namely muon 
capture of fi~ —>■ also becomes huge, and there is no chance of survival 
of captured muons left for the lepton number violating processes. Techni- 
cally, the relevant quantity in this case is the branching fraction of muon 
captured reactions instead of the absolute event rate, and for this quantity 



6 



the enhancement factor due to the capture into the atomic orbit cancels. 



We next discuss the processes off atomic electrons, ^ and (jH)). Both 
processes are quite unique in that the lepton number disappears in the final 
state. For the atomic target the threshold for both processes opens up at 
around ^ ^ SOGeV . The cross section for the process e + e is given by 



where /,r is the tt decay constant of order 90MeV, and s ~ 2me-E'e is the CMS 
energy squared. We have treated the atomic electron being at rest, since its 
momentum is much smaller than the incident beam energy E(.{~^ Tn^.). The 
important prefactor f 2^^"^ is numerically ~ 8x 10^^'' cm^ (|me^|/100eV^)^. 
The event rate assuming an effective flux of order 10^^ cm~^sec~^ and a heavy 
target of mass 500(7 and A ~ 2Z is 4 x 10^ {\me^\/lQQeVY /year . A non- 
trivial background of the usual electromagnetic origin is e~+n e^+j9+7r^ + 
TT^ + 7r+ with a missing p. The detector should be arranged hermetically not 
to miss particles produced in the final state, and the kinematical condition 
of the two-body process of ((Tj) should be maximally exploited to reject these 
backgrounds. 

The corresponding inclusive process, + ^ X , with X any hadronic 
state has a larger cross section of order, 10~^ G^lmeep , or 
^ 10-5Scm2(|mee|/100e1/)2(v^/100G'e1/)^ for ^ < 250GeV . A possi- 
ble experimental problem would be photons from 7r° decay which might be 
misidentified by electrons at high energies. Thus, it would be preferable, 
but difficult, to solely detect all charged modes not containing in the fi- 
nal state. Furthermore, for the atomic electron target a reasonable choice 
of the maximal ^/s would be of order IGeV, and even for this one needs a 
TeV electron LINAC. A more precise formula of the inclusive cross section 



For example, the branching fraction of the process jj, 
0[Gl\m^e?Z/Sn^R^] ~ 3 x 10-29|m^e/100eV|2[5]. 



e"*" is given by 




(13) 
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may be computed by using the simple quark model of 3 colors, to give at 
2mw > > 2m^ 



where -y/sxj is the invariant mass of the quark pair, Uidi. The integral of the 
invariant mass distribution /(xi , X2) is about 0.07. 

The neutrino beam is attractive from two reasons: the low background 
and a high intensity beam may become available in the proposed neutrino 
factory. The high intensity z/^ beam of energy of order SOGeV may also be 
of interest due to the possibility of the r-neutrino appearance experiment. 
A related tri-muon process to measure m^^ has also been considered in [0], 
which belongs to another class of effective operators, llUqq . Incidentally, 
we considered this class of operators along with yet another class of ////// , 
searching all promising processes of the lepton number violation. Due to 
the experimental difficulty of determining the sign of the lepton number of 
neutrinos the process of [9, , namely z/^+p — > fi^+fi'^+fi^+X , seems the only 
feasible process of both of these two classes for signatures of unambiguous 
lepton number violation. 

The cross section for (jSI) is just 1/4 times of the rate eq. lfT^ . neglecting 
the pion mass difference. A non-trivial background comes from the usual 
neutral current weak process, i/^ + n — z/^ + p + tt^ + tt'^ with a missing 
p. The true signal is kinematically selected by, for instance, plotting the 
invariant mass((p_ +Po)^) distribution of two pions (low energy 7r° may also 
help by requiring two photon showers of 7r° —>■ 77 decay) , which peaks around 
at \/2meEy for a given incident neutrino energy Ej^ . 

Finally, a low energy proton beam may be used in (jHl). The cross section 
at low energies for this process is O[10~™]cm^ for Z = 50. 

A meaningful limit of the absolute scale of neutrino masses is set by the 
current negative search of the neutrinoless double beta decay, which gives a 




(14) 
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limit PU], for the combination, 

\mee\ = \T.kU^k^k\ < 0.3eV . (15) 

This is a combination of the mass matrix element different from that appear- 
ing in the above process of (jD), namely mg^, but is identical to that 
of (0). For instance, with a phase cancellation the neutrinoless double beta 
decay and the related process ((Tj) might be suppressed, but the other pro- 
cesses involving might be relatively strong. On the other hand, WMAP 
has recently derived from a detailed fluctuation map of the microwave back- 
ground the following limit JT], ^21 'nT'max < 0.2eV . If this limit is to be 
taken at a face value, it may not be necessary to set a neutrino mass limit 
of order leV or larger. Nonetheless, we believe that terrestrial experiments 
are to be supplemented to exclude with certainty even this range of neutrino 
masses. Complementary and redundant information on neutrino mass ma- 
trix elements mai3 is crucial for further understanding of the lepton, hence 
the GUT sector. Thus, despite very small rates experimental search for these 
new reactions is very welcome. 

A precise relation between the parameters of the neutrino oscillation and 
the lepton number violating process is given by, 

E^m„^m^^ = T^kUakU^k ^l- (16) 

The quantity in the right hand side can be determined by a precision data of 
the spectral distortion of K2K and Kamland oscillation experiments, when 
it is collaborated with the absolute mass scale determined by a direct mass 
measurement. The (absolute values of) m^^ in the left hand side, on the other 
hand, could be measured through the lepton number violating processes. In 
the physical observable in Eq. ()16|) . however, the information of Majorana 
phases is lost, due to the product of mass matrix m and its hermitian conju- 
gate m) . This is why we should study the lepton number violating processes. 
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which are caused by the mass matrix m itself containing the information of 
the CP violating Majorana phases. 

To clarify this we write the unitary matrix U as 



where Umns is the Maki-Nakagawa-Sakata matrix with 1 CP violating phase 
and P = diag(l, e*°, e*^), with a and (3 being CP violating Majorana phases. 
We easily find that the observable in (fTB|) can be written as 
{Umns " diag{ml, ml, 171"^) ■ Uljp^g)ai3, where the matrix P disappears. On 
the other hand, from Eq.(3) we realize that ma/3 itself does depend on the 
Majorana phases: 



Having valuable information on Umns aiid the mass eigenvalues mi, m2, m^ 
from the neutrino oscillation experiments, by use of Eq. (jl8p we can evalu- 
ate each of {majsl and therefore the feasibility of lepton number violating 
processes. In the approximation of neglecting ^13 and taking the maximal 
mixing for the atmospheric neutrino oscillation, the unitary matrix Umns is 
expressed by 



where s© = siuOq and c© = cos 6'©. Substituting this expression for U m ns iii 
Eq. lfTHj) . we can readily compute \mai3\ ■ We generally find that = Im^^l 

and |me^| = |mer|. To get further results, we consider three typical cases of 
neutrino mass eigenvalues, i.e., the hierarchical (H), the inverted hierarchical 
(IH) and the quasi-degenerate (QD) cases. For the case H, utilizing m^ ~ 
\l^mlt^ » ma ~ Aml^i » mi, we find 



U — UmnsP, 



(17) 



ma/3 = {Umns ■ diag(mi, mse^'", mge^*'^) ■ Ul^^g)ap. 



(18) 



Umns 



/ Cq Sq o \ 

-Sq/V2 Cq/V2 -I/V2 



(19) 




10 



2V2 



2 V — til l 1 1 • 2 

1 / . O 1 O 



I^^MtI - |\/^"^atm - 24c2(a-/3)VAm2^; . (20) 

Therefore, the niee which is the effective mass for the neutrinoless double beta 
decay is small and also the contribution of the Majorana phases is suppressed. 



For the IH case, by using mi ~ 777-2 ~ y ^^atm >> we find 

mi^l- sIqsI > mi\c2Q\ ~ ^\/Amltm ~ 0.03eV , 



777p 



e/u| — /ri V itm 1 



a/2 

In this case, the effective mass of neutrinoless double beta decay is bounded 
from below and the decay will be measurable in high precision experiments, 
thus providing a useful information on the Majorana phase a. For the QD 
case, by using 7771 ~ 7772 — 7773 > Amatm: we find 

\mee\ ^ mi^l- sIqsI > 0.03eV , 

71" 



\m^^,\ ~ mi^l- (4c|s2 + sis} + c|s: 



^ mi^i- (44^^ + 44 + 44-/9)- (22) 



Therefore, again we have an enough chance to measure the effective mass 
of the neutrinoless double beta decay. In addition, from the experiments to 
measure the other matrix elements we may be able to obtain the information 
on both of two Majorana phases a and p. 

We thus learn that if the neutrinoless double beta decay unfortunately 
ends up with null result of precision ofO[0.03]eV, only the normal hierarchical 
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mass pattern is allowed. In this case the lepton number violation can be 
verified in terrestrial laboratories only by lepton number violating processes 
discussed in the present work. 

The interesting question is how CP violation, relevant for the leptogene- 
sis, may be examined by studying the CP violating phases in the Majorama 
mass matrix of left-handed neutrinos, which are obtainable from the neu- 
trino oscillation experiment, handled by Umns, and lepton number violating 
processes, as discussed above. To get the answer, we assume the seesaw mech- 
anism as the origin of the Majorana neutrino masses. In this mechanism the 
relevant quantity for the leptogenesis is m^m/j, with m^j being the Dirac 
mass matrix in the basis where the Majorana mass matrix of right-handed 
neutrinos is diagonalized. The Dirac mass matrix mD can be expressed as 
m£) = U\/DyR\/DR, where U is the mixing matrix given in Eq. lfTTj) . D^, 
and are the diagonal matrices whose eigenvalues are the observable small 
Majorana masses of left-handed neutrinos and the right-handed Majorana 
masses, respectively/,, i? is a complex "orthogonal" matrix, R = 1, with 3 
independent CP phases, which is otherwise arbitrary. It may be worthwhile 
to observe m^m^) = ^ DpiRDj^R'^ \/ Dr. That is, the mixing matrix U is not 
directly related to the leptogenesis. We realize that not only the phase in 
Umns, but also the Majorana phases disappear in m^mi). It, however, will 
be generally possible to relate the CP violating phases in U to those in R, 
responsible for the leptogenesis, once some complementary information, such 
as the absolute values of the elements of is obtained. It may also be 
possible that a relation between the phases of U and those of R is naturally 
realized in a sophisticated concrete model of neutrino mass generation. 

A road map towards experimental verification of leptogenesis may not 
be easy to draw. A complete determination of the Majorana mass matrix 
including the CP phase may require a variety of lepton number violating 
reactions such as nn eepp (neutrinoless double beta process), e~ — > /x"*", 
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Vfj^ — > iJi~ jjL^ jji^ ^ and pp — > T'^T'^nn. Even this much is not sufficient; one 
further needs a handle to the mass scale of the heavy Majorana particle in 
the seesaw mechanism to compare with leptogenesis calculation. 

We finally make a short comment on processes not discussed so far. The 
binding system of stable atoms (e~ + nucleus) has been considered since in 
this case a shorter time scale different from the reaction time is available 
(hke in the case of neutrinoless double beta decay), but we did not find a 
useful system due to the difficulty of overriding the barrier of nuclear binding. 
For instance, the He atom cannot spontaneously decay emitting e"*". As to 
the decay process, — > /i"*" + X , — > /i+ + yU+ + X etc. with X any 
hadronic state is conceivable, which has however a branching ratio of order, 
{GF\mij_r\^TY ~ 10"^^(|m^^|/100eV)^ for r, too small even for the 
tau-charm factory. (The current upper limit is of order 10~^.) 

Thus, to the best of our knowledge, the processes discussed in the present 
work are the best candidates to explore the Majorana nature and its strength 
of the neutrino mass matrix, beyond the neutrinoless double beta decay. If 
these tiny rates are experimentally falsified by larger rates, it definitely im- 
plies a new source of lepton number violation besides the Majorana neutrino 
mass. 
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